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Abstract. We compare four different methods to calculate radiative 
transfer through a magnetized stellar atmosphere, and apply them to the 
case of magnetic white dwarfs. All methods are numerically stable enough 
to allow determination of the magnetic field structure, but distinctions 
between faster, simplifying, methods, and elaborate, but more CPU-time 
consuming, methods, can be made. 

About 3% of all known white dwarfs have strong magnetic fields between 
10 6 and 10 9 Gauss (Wickr amasinghe & Ferrario 2000, Jordan 2001). The de- 
tailed surface structure of the magnetic field can be inferred from time-resolved 
spectro-polarimetric observations, which give integrals of the polarized radiation 
over the stellar surface visible at a given rotational phase. Because the magnetic 
field varies in strength and direction, theoretical calculations for polarized radia- 
tive transfer have to be performed at many (typically 1000) different points on 
the surface. In order to ensure numerical accuracy and efficiency, we have tested 
four different methods from the literature for solving the radiative transfer equa- 
tions in magnetized white dwarf atmospheres. Note that the ALI (Accelerated 
A-Iteration) method is described in more detail by Deetjen et al. in these pro- 
ceedings. 



1. The Radiative Transfer Equations 

The polarization properties of light are described by the four Stokes parameters 
/, Q, V, and U. I is the intensity. Linear polarization is described by Q and 
U, and circular polarization by V. In general, the absorption coefficients k\ 
and K r for left- and right-handed circularly polarized light, and k p for linearly 
polarized radiation traveling perpendicular to the magnetic field, are not equal at 
a given wavelength, m, k p , and K r correspond to transitions where the magnetic 
quantum number changes by Am = —1,0, and +1. The absorption coefficients, 
normalized to the Rosseland mean opacity (kr ss), are defined by r\ v = k p /kr oss , 
Vl = ^//^Ross ; and T] r = kv/kross- These absorption coefficients are combined 
into rji = i^sin 2 ^ + \{r]i + rj r )(l + cos 2 -0), VQ = Ivp ~ + Vr)sin 2 ip, 
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and i]v = \{f]r — r)i)cosip. Here ip denotes the angle between the magnetic 
field direction and the line of sight. When the magneto-optical parameters pr 
(Faraday rotation) and pw (Voigt effect) originating from spectral lines and free 
electrons in a magnetic field (see Jordan et al. 1992) are taken into account, the 
three radiative transport equations of Unno (1956) expand into four equations 
(Beckers 1969): 

VlTZ = Vi(I ~ B) + r] Q Q + W V, (1) 



dr 

^ = vq(i-b) + V iQ + prU, 

dll _ 
p— = p R Q + rijU - p w V, 
dr 

dV 

P— = r) V {I - B) + p w U + 7] V V. 
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The optical depth scale is dr = —Kn oss dz', and p = cos#, where is the 
angle between the normal to the surface and the line of sight (along which the z' 
axis of the local coordinate system on each surface element points). The x' axis 
is defined by the direction of the projection of the magnetic field vector onto 
the plane of the sky. yi completes the right-handed coordinate system. B is the 
source function (the Planck formula in the case of LTE) . 

In the case of large Faraday rotation, radiative transfer can be described 
in the framework of generalized Stokes parameters, with I + being the ordinary 
mode and i_ being the extraordinary mode. For the ordinary mode the electric 
field vector, the magnetic field, and the direction of propagation are in one plane; 
for the extraordinary mode the plane of propagation and the electric field vector 
is perpendicular to the magnetic field. The phase relation between the two 
modes is given by I c = ^/I+I^ cos 5, and I s = \/T+TZ sin 5, which correspond to 
two polarization ellipses, the so-called normal modes. 

We have used three different methods to solve the radiative transfer equa- 
tions formally: (a) The slightly modified Martin und Wickramasinghe (M&W, 
1979), (b) an accelerated lambda iteration scheme (ALI), (c) an approximation 
for large Faraday rotation (APPROX), or (d) an analytical procedure using 
matrix exponential function (MATEXP). 

1.1. Algorithm by M&W 

The solution is based on the assumption that the source function is linear in the 
optical depth and that between two successive depth points the Stokes parame- 
ters can be described by 
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(5) 



This assumption is inserted into the four radiative transfer equations, using an 
Unno (1956) condition at the inner boundary. Then the coefficients I a , . . . V c i, 
as well as dj, are evaluated by comparing the constant, linear, and exponential 
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terms on both sides of the equations. One obtains finally a relatively simple 
scheme, by which the transport equations can be solved from the inner depth 
points to the outside. 

1.2. ALI method 

This powerful and very flexible method is described in more detail by Deetjen 
et al. in these proceedings. 

1.3. APPROX method 

This method uses generalized Stokes parameters (Ramaty 1969). Of particular 
interest is the limiting case of large Faraday rotation, which is always valid in 
the case of magnetic white dwarfs. This means that the rotation between two 
depth points is larger than 2tt, and the phase relation between the polarization 
ellipses I s and I c can assume arbitrary values; I s = I c = on the average. Then 
the transfer equations decouple (Ramaty 1969): 

"IF = a± (I " 7± ) ' (6) 

According to Vath & Chanmugam (1995) the formal solution (using a linear 
approximation of the source function) of this equation can be written as 

T i+1 e' a + AT ^) - (1 - e -«+Ar/M)" 



r i+1 e- a - AT ^) - (1 - e - a - AT/ v) 



+/ + (T, ±1 )e- a + A ^ 

i-(n) = \a B (i-e- a ~^/n\—U + {n 

I 2 a_ l 

+I_(r t±1 )e- a - AT ^ 



1.4. MATEXP method 

A direct analytical way to solve the radiative transfer equations uses matrix 
exponential functions. After Dittmann (1995) the radiative transfer equation 

dl/dr = r)(I — S), where / = (I, Q, V, U), can be written as 



1(0) = exp(— ijt) 



I(t) — T) I exp(—rjs)S(s)ds 



(7) 



where the matrix exponential function is defined as = X)«^=o r 7™/ n " The 

3 

infinite series can be rewritten as e^ T = y~2 q(t)M ! . The matrix M is given by 

i=0 

M = T] — r//E. We assume a linear source function between successive depth 
points: S = (a + 0t, 0, 0, 0). If we define x = \\J%Jli + (%), y = |y / 2(^|^p|), 

z = -iy, pi = \Jvi + 4 (VQPw + WPr) 2 , V5 = \[tf Q + Vv ~ Pr~ Pw and take 
into account that y can be become imaginary if magneto-optical parameters are 
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present (a case that has been neglected by Dittmann 1995), then the coefficients 
are given by: 



co 



Cl 



C2 



C3 



e Vir 
e Vir 
e Vir 

7f 



x 2 cos(zr) — y 2 cosh(xr) 



x 2 y 2 

— sin(zT) sinh(xr) 

Z X 



[cosh(xr) — cos(zr)] 

— sinh(xr) sin(zr) 

x z 



(8) 
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For large magneto-optical parameters lim^o sinh(xr)/a; = r, and the expres- 
sions can be somewhat simplified. 



2. Conclusion 

In all tests, the resulting spectra and polarization turned out to be very similar 
for the four different approaches. Small numerical instabilities in the case of 
the M&W and the MATEXP algorithms lead to small oscillations of the Stokes 
parameters (particularly Q and V) in the vicinity of spectral lines; thsese are 
caused by the terms containing sin, cos, sinh, or cosh. The ALI method, having 
the potential to account for NLTE and scattering is numerically very stable. 
However, the APPROX method is about a factor of 11 faster and very stable and 
should therefore be preferred for practical purposes. For final fits to observations 
it is always possible to use the ALI method for comparison. 
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